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Accurate infrared (IR) spectroscopy of protonated water clusters prepared in 
the gas phase has become possible in recent years, [IJGJ El HI El E] opening the 
door to a deeper understanding of the properties of aqueous systems and the 
hydrated proton, which are of main interest in central areas of chemistry and 
biology. Several computational studies have appeared in parallel, providing 
a necessary theoretical basis for the assignment and understanding of the 
different spectral features. El El El HO] It has been recently demonstrated 
that the HsO^" motif, also referred to as Zundel cation, plays an important 
role in protonated water clusters of 6 or more water molecules and as a 
limiting structure, together with the Eigen cation (H9O4"), of the hydrated 
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proton in bulk water. [U El [TT] The importance of the hydrated proton and 
the ammount of work devoted to the problem contrast with the fact that 
the smallest system in which a proton is shared between water molecules, 
H5O2", is not yet completely understood, and an explanation of the most 
important spectral signatures and the associated dynamics of the cluster is 
lacking. 

In this letter we report the simulation of the IR linear absorption spec- 
trum of the cation in the range between to 4000 cm -1 by state- 
of-the-art quantum-dynamical methods, we discuss the spectral signatures 
in terms of the underlying couplings and dynamics of the different degrees 
of freedom and compare our results to recent, accurate experiments on this 
system. For the first time the doublet-peak feature at around 1000 cm -1 
is fully reproduced, analyzed and assigned. The doublet is found to arise 
from the coupling between the proton-transfer mode and low frequency, 
large-amplitude displacements of both water molecules. Predictions are also 
made for the lowest frequency part of the spectrum, which has not yet ex- 
perimentally been accessed. Several important features of the system are 
analyzed for the first time, namely the degeneracy of some of the vibrational 
levels and the extreme anharmonicity of the wagging motions and relative 
internal rotation of the two water molecules. In doing so, we do not resort 
to any low-dimensional model of the system, but we treat it in its full di- 
mensionality, i. e., 3iV — 6 = 15 active internal coordinates (15D). The use 
of full dimensionality is found crucial in the reproduction of the complete 
absorption spectrum and dynamics. Our study provides a picture of the 
HsO^ system, extendable to larger aggregates, in which the clusters have to 
be viewed as highly anharmonic, flexible, multi-minima, coupled systems. 
From a methods perspective, we show that a full quantum-dynamical de- 
scription of such a complex molecular system can still be achieved, providing 
explicative and predictive power and a very good agreement to available ex- 
perimental data. In this respect, the reported simulations set a new state of 
the art in quantum dynamically describing an anharmonic, highly coupled 
molecular system of the size of the HsO^ cation. To account for the inter- 
atomic potential and the interaction with the radiation we make use of the 
potential energy surface (PES) and dipole-moment surfaces recently devel- 
oped by Bowman and collaborators, [9] which constitute the most accurate 
ab initio surfaces available to date for this system. 

The IR predissociation spectrum of the HsO^ cation has been recently 
measured in argon-solvate [3] and neon- and argon-solvate [6] conditions. 
It is expected that the photodissociation spectrum of the HsO^-Nei com- 
plex is close to the linear absorption spectrum of the bare cation. [5] This 
spectrum features a doublet structure in the region of 1000 cm -1 made 
of two well-defined absorptions at 928 cm -1 and 1047 cm -1 . This dou- 
blet structure was not fully understood, although the highest-energy com- 
ponent was assigned to the asymmetric proton-stretch fundamental based 
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on the quantal calculations on the IR spectrum of H5O2 of Bowman and 
coworkers. [6 J Low-frequency modes may also play an important role in com- 
bination with the proton-transfer fundamental. Such a possibility has been 
already suggested, [21 [HJ [TO] but just which modes would participate in such 
combinations, and how, is still a matter of discussion. 

The Hamiltonian used in the simulation of HsO^ is expressed in a set 
of polyspherical coordinates based on the Jacobi vectors in Figure l.|12j It 
is found that only after the introduction of such a curvilinear set of coordi- 
nates an adequate treatment of the anharmonic large-amplitude vibrations 
and torsions of the molecule becomes possible. The kinetic energy operator 
is exact for J = 0, and the derivation of its lengthy formula (674 terms) 
will be discussed in a forthcoming publication. The correctness of the op- 
erator implemented was checked by comparison with data generated by the 
TNUM program. [13J The internal coordinates used are: the distance be- 
tween the centers of mass of both water molecules (R), the position of the 
central proton with respect to the center of mass of the water dimer (x,y,z), 
the Euler angles defining the relative orientation between the two water 
molecules (wagging or pyramidalization: 7 a ,7f>; rocking: (3 a , fib', internal 
relative rotation: a) and the Jacobi coordinates which account for the par- 
ticular configuration of each water molecule (ru^fy, r 2 ( a! &), #( a ,6))) where r\ x 
is the distance between the oxygen atom and the center of mass of the corre- 
sponding H2 fragment, r^ x is the H-H distance and X is the angle between 
these two vectors. These internal coordinates are body-fixed (BF) ones, 
where the water-water distance vector R points along the BF z-axis. These 
coordinates have the great advantage of leading to a much more decoupled 
representation of the PES than a normal-mode based Hamiltonian. The 
quantum-dynamical problem is solved in the time-dependent picture using 
the multiconfiguration time dependent Hartree method (MCTDH). [14l [T5] 
The potential energy surface has been represented by a cut high-dimensional 
model-representation (cut-HDMR).[TO| I17j 

In Figure 2 probability-density projections on the wagging coordinates 
are shown for the ground vibrational state (go), as well as for one of the two 
fundamental states (w\ a ,Wxb) of the wagging modes, which are degenerate 
vibrational states with an energy of 106 cm -1 . State w% is shown in Figure 
2c and it will play a major role due to its coupling to the proton-transfer 
mode, as will be discussed later. 

The probability-density of the wagging coordinates in go (Figure 2a) 
presents four maxima in which the wagging angle is about 30 degrees with 
respect to the planar conformation for each water molecule. The probabil- 
ity for one or both water molecules to be found in a planar conformation 
is almost as high as the probability to be found pyramidal. This means 
that HsO^ interconverts already at T = 0, due to zero-point energy, be- 
tween equivalent absolute minimum-energy structures in which both water 
molecules are found in a pyramidal conformation. Four equivalent minimum- 
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energy structures are accessible through wagging motions. The number of 
accessible equivalent minima at T = doubles to eight since the relative 
rotation of both water molecules (a coordinate) has also been found to be 
allowed through a low energy barrier. 

The energies of the next three wagging-mode states (u>2, 1^3,^4) are, re- 
spectively, 232, 374 and 422 cm -1 . In a harmonic limit these states can be 
represented by kets |11), (|20) - |02))/ v / 2and (|20) + |02))/ v / 2, respectively, 
where the | ab) notation signifies the quanta of excitation in the wagging mo- 
tions of waters a and b. The degeneracy between W2, u>3 and is broken 
due to anharmonicity. In harmonic approximation the energies of the two 
lowest wagging- fundamentals w\ a and wib are about 300 cm -1 larger than 
our result and do not account for their degeneracy, since harmonic modes are 
constructed taking as a reference only one of the equivalent absolute minima. 
The system, however, interconverts between eight equivalent Ci structures 
and other stationary points through low-energy barriers (wagging motions 
and internal rotation) , which leads to a highly symmetric ground-state wave- 
function. Other vibrational states have been computed which are related to 
the internal rotation, rockings and water-water stretching modes. They will 
be reported and discussed in a forthcoming publication. 

Figure 3 presents the IR predissociation spectrum of the HsO^-Ne com- 
plex [6] and the simulated spectrum of HsO^ in the range 700-1900 cm -1 . 
The simulated spectrum is obtained in the time-dependent picture by Fourier 
transformation of the autocorrelation of a dipole-operated intial state: [T8] 

Eq. 1 HERE 

where Eq is the ground-state energy and |^,o) = Al^o)- The simulated 
spectrum shows a good agreement with the experimental spectrum. The 
agreement on the doublet structure around 1000 cm^ 1 is very good, and 
the position of the doublet at 1700 - 1800 cm -1 is also in good agreement, 
despite the relative intensities being larger in the simulation. 

The simulated spectrum in the range between and 4000 cm -1 is de- 
picted in Figure 4. The region below 700 cm -1 has not yet been accessed 
experimentally. Direct absorption of the wagging motions, excited by the 
perpendicular components of the field, appears in the range between 100 - 
200 cm -1 . The doublet starting at 1700 cm -1 is clearly related to bending 
motions of the water molecules, but its exact nature is still to be addressed. 
The simulated spectrum also shows the absorptions of the OH stretchings 
starting at 3600 cm -1 . 

The doublet absorption at around 1000 cm -1 and the related underlying 
dynamics deserve a deeper analysis. Due to the high density of states, it 
was not possible to obtain the fully converged states, but reasonably good 
approximations to the wavefunctions of the low-energy (l^), 930 cm -1 ) 
and high energy (|\^), 1021 cm -1 ) eigenstates of the doublet were com- 
puted. Even though these wavefunctions contain all the possible informa- 
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tion on the two states, their direct analysis becomes complex due to the 
high dimensionality of such objects. In order to obtain a fundamental un- 
derstanding of the observed bands, zeroth-order states and \$r W3 ) 
were constructed, where |$ z ) is characterized by one quantum of excitation 
in the proton-transfer coordinate whereas \&r >W3 ) by one quantum in the 
water-water stretch and two quanta in the wagging motion. They were con- 
structed by operating with z on the ground state: |<3? 2 ) = z\^o)N, where N 
is a normalization constant, and by operating with (R — Rq) on the third 
excited wagging state W3: |$r jU , 3 ) = (R — Ro)\^ W3 )N, respectively. The 
two eigenstates corresponding to the doublet were then projected onto these 
zeroth-order states. The corresponding overlaps read: |( ( I ) 2 | l I'^)| 2 = 0.20, 
K^l^)! 2 = 0.53 and |<**|*5>| 2 = 0.48, |($^ 3 |^)| 2 = 0.12. One 
should take into account that these numbers depend on the exact defini- 
tion of the zeroth-order states, which is not unique. However, they provide 
a clear picture of the nature of the doublet: the low-energy band has the 
largest contribution from the combination of the water-water stretch and 
the third excited wagging (see Figure 2c), whereas the second largest is 
the proton-transfer motion. For the high-energy band the importance of 
these two contributions is reversed. Thus, the doublet may be regarded as 
a Fermi resonance between two zero-order states which are characterized by 
(R, W3) and (z) excitations, respectively. The reason why the third wag- 
ging excitation plays an important role in the proton-transfer doublet is 
understood by inspecting Figure 2c and Figure 5. The probability density 
of state 1113 has four maxima, each of which corresponds to a planar con- 
formation of H2O-H 4 " (H30 + character) for one of the waters, and a bend 
conformation (H2O character) where a lone-pair H2O orbital forms a hydro- 
gen bond with the central proton. When the proton oscillates between the 
two waters, the two conformations exchange their characters accordingly. 
Thus, the asymmetric wagging mode (11)3, 374 cm -1 ) combines with the 
water- water stretch motion (R, 550 cm -1 ) to reach an energy close to the 
natural absorption-frequency of the proton transfer. As a consequence, the 
low-frequency wagging (or pyramidalization) motion of the water molecules 
becomes strongly coupled to the higher frequency, spectroscopically active 
proton-transfer motion and this coupling leads to the characteristic doublet 
feature of the IR spectrum. 

In conclusion, we report a simulation of the dynamics and IR absorption 
spectrum of the HsO^ cation by quantum-dynamical methodology in the full 
spectral range 0-4000 cm" 1 . The spectrum is directly comparable to avail- 
able and future experiments on this system. We discuss some important fea- 
tures of the protonated water dimer which have remained until now elusive 
due to strongly anharmonic, large-amplitude motions. The floppy, multimin- 
ima nature of the cluster is presented and analyzed and the doublet-band 
absorption around 1000 cm" 1 is fully reproduced and explained in terms 
of coupling of the proton-transfer motion to wagging torsions of the water 
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moieties. These calculations constitute an avenue for a detailed quantum- 
dynamical description of larger clusters and provide important fundamental 
information on the spectroscopy and dynamics of protonated aqueous sys- 
tems and the hydrated proton. 
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Figure Captions: 



Figure 1: Set of Jacobi vectors in terms of which the kinetic energy of the 
system is expressed. The set of internal coordinates used corresponds to the 
length of these vectors and relative angles. The z direction of the central 
proton is parallel to R. 

Figure 2: a) Probability density of the ground vibrational state, b) first 
and c) third wagging-mode states projected onto the wagging coordinates 
7 a and j b . 

Figure 3: a) Predissociation spectrum of the HsO^-Ne complex;[6] b) 
quantum-dynamical simulation. The resolution at which different peaks 
of the spectrum are resolved is given by the Fourier Transform. Due to the 
finite propagation time a finite resolution of about 30 cm -1 is obtained. 
Figure 4: Quantum-dynamics simulated spectrum in the range between 
and 4000 cm -1 . Absorption is given in absolute scale in mega-barns (Mb). 
Figure 5: Two most important coupled motions responsible for the doublet 
peak at 1000 cm" 1 . 
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